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We report the study of the helicity driven photocurrents in graphene excited by mid-infrared light 
of a CC>2-laser. Illuminating an unbiased monolayer sheet of graphene with circularly polarized 
radiation generates - under oblique incidence - an electric current perpendicular to the plane of 
incidence, whose sign is reversed by switching the radiation helicity. We show that the current is 
caused by the interplay of the circular ac Hall effect and the circular photogalvanic effect. Studying 
the frequency dependence of the current in graphene layers grown on the SiC substrate we observe 
that the current exhibits a resonance at frequencies matching the longitudinal optical phonon in 
SiC. 

PACS numbers: 73.50.Pz, 72.80.Vp, 81.05.ue, 78.67.Wj 



I. INTRODUCTION 

Recently graphene has attracted enormous attention 
because its unusual electronic properties make possi- 
ble relativistic experiments in a solid state environment 
and may lead to a large variety of novel electronic de- 
vices [2h3- One of the most interesting physical aspects 
of graphene is that its low-energy excitations are mass- 
less, chiral Dirac fermions. The chirality of electrons in 
graphene leads to a peculiar modification of the quan- 
tum Hall effect 0, [?J , and plays a role inphase-coherent 
phenomena such as weak localization [1, Q . Most of cur- 
rent research in this novel material are focused on the 
transport and optical phenomena. In our recent work, 
we reported on the observation of the circular ac Hall 
effect (CacHE) [To[ which brings the transport and opti- 
cal properties of graphene together: In CacHE an electric 
current, whose sign is reversed by switching the radiation 
helicity, is caused by the crossed electric and magnetic 
fields of terahertz (THz) radiation. The photocurrent is 
proportional to the light wavevector and may, therefore, 
also be classified as photon drag effect [10l - [l6l | . Classical 
theory of CacHE, well describing the experiment at THz 
frequencies, predicts that for lot 3> 1, with lo being the 
radiation angular frequency and r momentum relaxation 
time of electrons, the ac Hall effect is suppressed. 

Here we demonstrate, however, that helicity driven 
photocurrents can be detected applying a mid-infrared 
CO2 laser operating at much higher light frequencies 
where the condition lot ^> 1 is satisfied. Our results 
show that in this case, due to the fact that the classi- 
cal CacHE is substantially diminished, much finer effects, 
such as circular photogalvanic effect (CPGE), well known 
for noncentrosymmetric bulk and low dimensional semi- 
conductors [TBI [l7l - r20j , become measurable. We present 
a phenomcnological and microscopic theory of photocur- 



rents in graphene and show that the experimental proof 
of the interplay of CacHE and circular PGE of compa- 
rable strength comes from the spectral behavior of the 
photocurrent. 

Our experiments demonstrate that variation of the ra- 
diation frequency may result in an inversion of the pho- 
tocurrent sign. We show that the light frequency, at 
which the inversion takes place, changes from sample to 
sample. Tuning the radiation frequency in the operation 
range of a mid-infrared CO2 laser we also observed a 
resonant-like behaviour of the photocurrent in graphene 
grown on the Si-terminated face of a 4H-SiC(0001) sub- 
strate: its amplitude drastically increases at frequency 
/ = 29.2 THz (A = 10.26 fim). The microscopic origin of 
the resonant photocurrent is unclear, but we show that 
its position is correlated with the high frequency edge of 
the reststrahlen band and, correspondingly, to the energy 
of the LO phonon in 4H-SiC. Besides the helicity driven 
electric currents we also present a detailed study of a pho- 
tocurrents excited by unpolarizcd and linearly polarized 
light, also observed in our experiments, and discuss their 
origin. 



II. EXPERIMENT 

The experiments were carried out on large area 
graphene monolayers prepared by high temperature Si 
sublimation of semi-insulating silicon carbide (SiC) sub- 
strates [2l[. The samples have been grown on the Si- 
terminated face of a 4H-SiC(0001) substrate. The reac- 
tion kinetics on the Si-face is slower than on the C-face 
because of the higher surface energy, which helps ho- 
mogeneous and well controlled graphene formation [22| . 
Graphene was grown at 2000°C and 1 atm Ar gas pres- 
sure resulting in monolayers of graphene atomically uni- 
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FIG. 1: Photocurrent j normalized by the light intensity I as 
a function of the angle <p defining radiation helicity. Here 
j(ip) is measured at room temperature applying radiation 
with hui = 133.4 meV (A = 9.27 fim). Open and full circles 
show the longitudinal, j x , and transverse, j y , photocurrents 
measured at oblique incidence (8o = —30°) along and per- 
pendicular to the light propagation, respectively. Triangles 
demonstrate that the photoresponce vanishes at normal inci- 
dence (#o = 0°). Lines show fits according to Eqs. fTJ, {2} 
obtained using only photocurrent magnitudes as fitting pa- 
rameters [see also Eqs. ((3|, Q and ((H), ©]. The inset shows 
the experimental geometry, the plane of incidence of the ra- 
diation and the arrangement of contacts (black dots) at the 
edges of graphene. The ellipses on top illustrate the polariza- 
tion states for various tp for light incident on the sample as 
seen along the propagation direction. 
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FIG. 2: Angle of incidence dependence of the various 
photocurrent contributions detected in the transverse (up- 
per panel) and longitudinal (lower panel) geometries. Here 
j y ,A(9o), 3v,b{0o), jx,B'(8o), and j x ,c(&o) are obtained by 
measuring the helicity dependence of the photocurrent and 
fitting it by the Eqs. (O and © [see also Eqs. ©, @ and 
©]. The solid lines are fits after j oc 6q. The inset shows 
the experiment geometry. 



form over more than 1000 /im 2 , as shown by low-energy 
electron microscopy [23[. Four contacts have been cen- 
tered along opposite edges of 5 x 5 mm 2 square shaped 
samples by deposition of 3nm of Ti and 100 nm of Au 
(see inset in Fig.[lJ. The measured resistance was about 
2 kfi. From low-field Hall measurements, the manufac- 
tured material is n-doped due to the charge transfer from 
SiC 1221 24[ . We used two layers of non-conductive poly- 
mers [25j to protect graphene samples from the undesired 
doping in the ambient atmosphere and to control carrier 
concentration in the range (3 to 7)xl0 12 cm" 2 , mobility 
is of the order of 1000cm 2 /Vs and the Fermi energies 
Ep ~ 300 meV. All parameters are given for room tem- 
perature. 

To generate photocurrents we applied mid-infrared ra- 
diation of tunable C02-lasers with operating spectral 
range from 9.2 to 10.8 ^im (32.6 THz< / < 27.8 THz) 
correspond ing to photon energies ranging from 114 to 
135 meV [la]. For these wavelengths the conditions 
hui < Ef and lot ^> 1 hold. Two laser systems were used; 
a medium power Q-switched laser with the pulse dura- 
tion of 250 ns (repetition frequency of 160 Hz) and low 
power continuous- wave (cw) laser modulated at 120 Hz. 
The samples were illuminated at oblique incidence with 
peak power, P, of about 500 W and about 0.1 W for Q- 



switched and cw laser, respectively. The radiation power 
was controlled by photon drag detector (26| and/or MCT 
detector. The radiation was focused in a spot of 1 mm 
diameter being much smaller than the sample size even 
at oblique incidence [27j . The initial laser radiation po- 
larization vector was oriented along the a;-axis. Applying 
Fresnel A/4 rhomb we modified the laser light polariza- 
tion from linear to elliptical. The helicity P C i rc of the 
light at the Fresnel rhomb output was varied from -1 (left 
handed circular, er_) to +1 (right handed circular, er + ) 
according to P C irc = sin 2tp, where tp is the azimuth of 
Fresnels romb. Angle tp = corresponds to the position 
of the Fresnel rhomb when its symmetry plane is oriented 
perpendicular to the y-axis. The polarization ellipses for 
some angles <p are shown on top of Fig. [TJ 

The geometry of the experiment is sketched in the in- 
set in Fig. [T] The incidence angle #0 w &s varied between 
—30° and +30°. In our experiments we used both trans- 
verse and longitudinal arrangements in which photore- 
sponse was probed in directions perpendicular and paral- 
lel to the light incidence plane, respectively (see insets in 
Fig. [T] and (2). The photosignal is measured and recorded 
with lock-in technique or with storage oscilloscope. The 
experiments were carried in the temperature range from 
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FIG. 3: Helicity driven photocurrents given by the coefficient 
A = j y ,A/(I9o) as function of c<jt. Solid curves show calcula- 
tions of the ac Hall effect. The results of calculations and the 
low frequency data (ojt < 0.4) are given after [!(}, see also 
Eqs. ©. 

4.2 K to 300 K. 

The signal in unbiased samples is observed under 
oblique incidence for both transversal and longitudinal 
geometries, where the current is measured in the direc- 
tion perpendicular and parallel to the plane of incidence, 
respectively. Figure Q] shows the photocurrent as a func- 
tion of the angle (p for these geometries. The current be- 
haviour upon variation of radiation ellipticity is different 
when measured normal and along to the light incidence 
plane. 

The photocurrent for the transversal geometry, j y , (see 
full circles in Fig. [IJ is dominated by the contribution 
proportional to the photon helicity P C irc — sin 2^; it re- 
verses when the light polarization switches from the left- 
handed (cp = 45°) to the right-handed (ip = 135°) light. 
The overall dependence of j y on ip is more complex and 
well described by 

3v =j y ,Asin-2(p + jy tB sm4(p + £, (1) 

where j y .A = AIOq and j y ,B = BIOq are the magni- 
tudes of the circular and linear contributions, respec- 
tively. Here I is the light intensity. It is noteworthy 
that the offset £ is detected only in some measurements; 
It is almost zero and is neglected in the analysis below. 
The fit to the above equation is shown in Fig. [1] by solid 
line. We emphasize, that exactly the same functional be- 
haviour is obtained from a phenomcnological picture and 
microscopic models outlined below. Note that for cir- 
cularly polarized light, the current is solely determined 
by the first term in Eq. ([T]), because the degree of linear 
polarization is zero and, in this case, the second term 
vanishes. Our experiments show that jy,A and j y .B are 
odd functions of the incidence angle 9q ; a variation of 9q 
in the plane of incidence changes the sign of the currents, 
which vanish for normal incidence, #o~0 (sec triangles in 
Fig. [TJ. This behaviour is illustrated by Fig. [2] showing 
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FIG. 4: Spectral dependence of j y .A obtained for circu- 
larly polarized light (tp = 45°) and two angles of incidence 
6o = ±30°. Upper and lower panels show the data for sam- 
ples I and 2, respectively. The inset shows the temperature 
dependence of j Vt A measured in sample 2. 



the angle of incidence dependence of the photocurrents 
Ja, Jb and jc determining the magnitudes of the circular 
photocurrent and that depending on the degree of linear 
polarization, respectively. 

In the longitudinal geometry (open circles in Fig. [!}, 
the current sign and magnitude are the same for left- 
handed to right-handed circular polarized light and its 
overall dependence on tp can be well fitted by 

jx = jx,B> COS Aip + j XiC , (2) 

where j x ,B' = B' 10$ and j x ,c = CI8q are the magnitudes 
of the linear and polarization-independent contributions, 
respectively. The fit after this equation is shown in Fig. Q] 
by dashed line. Like in transversal geometry the pho- 
tocurrent angular dependence is in agreement with the 
theory discussed below. 

Figure [3] shows spectral behaviour of the circular pho- 
tocurrent given by the coefficient A = j y ,A/I9o- I n this 
figure A is plotted as a function of lot for both graphene 
samples. Besides the data obtained for light with the 
photon energy exceeding 110 meV we included here the 
results obtained in the same samples but at much lower 
THz frequencies / < 4 THz with fko < 16 meV. The 
latter data as well as the calculated dependences of the 
ac Hall effect are taken from our previous work [10| . It 
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FIG. 5: Spectral behaviour of A = j Vl A/(I&o) (left panel) 
and C = jx,c/(I0o) (right panel) in the vicinity of reso- 
nance. Solid curves show the reflection of the sample. We 
note that measurements of the reflection from the graphene 
and the back side of the sample yield the almost the same 
result. The inset shows the reflection spectrum in a larger 
frequency range. 

is seen that in the second sample the theory of the ac 
Hall effect describes well the experiment in the whole fre- 
quency range, including the high frequency data. While 
the sign and the magnitude of the current in the first 
sample measured at low frequency edge of the CCVlaser 
operation also fits well to the smooth curve of the ac Hall 
effect (see open circles in Fig. [3]) at high frequencies we 
observed that the signal abruptly changes its sign with 
rising frequency. The observed spectral inversion of the 
photocurrent's sign reveals that only ac Hall effect can 
not describe the experiment. 

Figure 2] shows the results of the more detailed study of 
the circular photocurrent's frequency dependence. The 
data were obtained by using the whole accessible, but 
very narrow, operating range of the C02-laser (114 meV 
< hui < 135 meV). Full and open circles in this figure 
correspond to the data obtained for two opposite angles 
of incidence 8o = ±20°. It is seen that the detected in 
sample 1 reversal of the current direction takes place at 
Huiinv ~ 119 meV. Here, uii nv indicates the frequency 
of the sign inversion. The drastic difference in the pho- 
tocurrent's spectral behaviour detected for samples with 
almost the same mobility and carrier density but pre- 
pared not in the same growth circle we attribute to the 
change of coupling between graphene layer and the sub- 
strate. In fact, this parameter is crucial for the mecha- 
nisms of the photocurrent generation. It may be different 
from sample to sample and it is difficult to control. 

Besides the spectral inversion, we observe another re- 
markable feature of the photocurrent: in both samples 
we detected a resonance increase of the current mag- 
nitude at hot ~ 121 meV (see Fig. 2] and left panel 
in Fig. [5]) . Similar resonance-like behaviour is detected 
for the polarization-independent longitudinal photocur- 



rent (sec right panel in Fig. [5]). The position of the 
resonance corresponds to the longitudinal optical (LO) 
phonon energy in 4H-SiC. In order to prove this we mea- 
sured the sample reflection for the graphene and the sub- 
strate sides. The results for both sides almost coincide 
with each other: the reflection shows the reststrahlen 
band behaviour (see the inset in Fig. [3]) . Solid curves in 
the left and right panels in Fig. [5] show that the high fre- 
quency edge of the reststrahlen band, which corresponds 
to the LO phonon energy in 4H-SiC, coincide with the 
resonance position. The detailed study of the resonance 
photocurrent and its power dependence is beyond the 
scope of the present work. 

To summarize the experimental part we demonstrate 
that illumination of graphene monolayers by mid- infrared 
radiation at oblique incidence results in the generation of 
photocurrents. Their directions and magnitudes are de- 
termined by the polarization of the radiation. At the 
frequencies about / = 29.2 THz (hui = 121 meV) we ob- 
served resonance feature and sign inversion of the pho- 
tocurrent. The latter property is sample dependent. 

III. THEORY 

Below we present phenomenological analysis of the 
photocurrents in graphene as well as their microscopic 
models. We demonstrate that the experimentally ob- 
served incidence angle, linear polarization and helicity 
dependences of the photocurrents correspond to phe- 
nomenological models. The magnitudes of the photocur- 
rents and their polarization dependencies are also in good 
agreement with theoretical predictions. 

A. Phenomenological analysis 

The ideal honeycomb lattice of graphene is described 
by the point group Dq^ containing the spatial inversion. 
As a result, photocurrent generation is possible provided 
that the joint action of electric, E, and magnetic, B, 
fields of the radiation is taken into account or provided 
that the allowance for the radiation wave vector, q, trans- 
fer to electron ensemble is made. In the former case the 
Cartesian components of the current are proportional to 
the bi- linear combinations E a B*p, while in the latter case 
to the combinations q a EpE*. Here Greek subscripts enu- 
merate Cartesian components. For the plane wave its 
wave vector, electric and magnetic fields are interrelated, 
therefore, for the purposes of the phenomenological anal- 
ysis it is enough to express the photocurrent density via 
the combinations q a EpE* as [16[ 

jx/I = r lg > |2 + K|2 + T 2 q J e ^-^\ (3a) 

g g* _|_ g* g 

j y /I = T 2 q x v - - V - fiq x P c i Ic e z . (3b) 
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where x and y are the axes in the graphenc plane, and 
z is the structure normal, the radiation is assumed to be 
incident in (xz) plane, e is the unit vector in light propa- 
gation direction and e is the (complex) polarization vec- 
tor of radiation, P c i rc is the circular polarization degree 
and q is the radiation wave vector. Additional contri- 
butions to the photocurrents, involving z component of 
electric field are analyzed in Ref. [l6j]. These effects arc 
expected to be strongly suppressed in ideal samples and 
for moderate radiation frequencies. Expressions (J3|) can 
be rewritten via incidence angle, 9 0l and angle tp deter- 
mining the radiation helicity as Eqs. ((J), ©• It allows 
one to establish a link between phcnomcnological con- 
stants Ti, T2 and T\ and fitting parameters A, B and C 
used to describe the experimental data, see Figs. Q] and 
[21 Namely, at small incidence angles 

AocTi, BocT 2 , and C ozT x . (4) 

It follows from Eqs. ^ that photocurrent contains, 
in general, three contributions illustrated in Fig. |6l 
panels (a)-(c). First one, schematically illustrated in 
Fig. inja) results in the polarization-independent pho- 
tocurrent flowing along the light incidence plane. 

In accordance with the general line of the paper we pay 
special attention to the photocurrent contribution pre- 
sented in Fig.lHJb) where the generation of the transversal 
to the light incidence plane current is shown. This cur- 
rent component is dependent on the radiation helicity: by 
changing photon from right- to left- circularly polarized, 
current changes its direction. This is nothing but the 
CacHE uncovered recently in graphene Ref. In ad- 
dition, transversal photoresponse contains a component, 
being sensitive to the linear polarization of radiation, see 
Fig. He). 

The photocurrent components described by Eqs. ([3]) 
can also be qualified as photon drag effects [HI, HH since 
in their phenomenological description the photon wave 
vector is involved. The direction of the photocurrent 
changes its sign upon reversal of the incidence angle. The 
contributions given by Eq. (JHaJ) and the first term in the 
Eq. (|3b[) can be easily understood as transfer of linear 
momenta of photons to the electron system [2^] and is 
recently discussed for graphene [TfJ [3(|. The circular 
photon drag current described by the second term on the 
right hand side of Eq. (|3b[) is due to transfer of both lin- 
ear and angular momenta of photons to free carriers. The 
circular photon drag effect was discussed phenomenolog- 
ically jU HI] and observed in GaAs quantum wells in the 
mid-infrared range (33j and in metallic photonic crystal 
slabs (34| . We note, that while the microscopic descrip- 
tion of the circular photocurrent in graphene in terms 
of ac Hall effect is relevant to the relatively low radia- 
tion frequencies range at high frequencies all photocur- 
rent contributions can be conveniently treated in terms 
of photon drag effect. 

The real structures, however, are deposited on a sub- 
strate, which removes the equivalence of the z and — z 
directions and reduces the symmetry to the Cq v point 




FIG. 6: Schematic illustration of the possible contributions 
to the photon drag and photogalvanic effects. Panels (a)- 
(c): polarization independent, circular and linear photon drag 
effects, Eqs. ((3|. Panel (d) shows the experimental geometry. 
Panels (e)-(f): photogalvanic effects allowed by symmetry in 
graphene samples deposited on substrates. 

group. Such symmetry reduction makes photogalvanic 
effects possible. The photogalvanic effects give rise to 
the linear and circular photocurrents fill ]: 

j x /I = Xi exet+ 2 e>Z , (5a) 

■ / T + e V eZ I D - fKU\ 

Jy/ 1 = Xi 2 h XcPch c e x , (5b) 

described by two independent parameters xi an d Xc- 
Schematically, these contributions to the photocurrent 
are shown in Fig. [SJe), (f). It follows from Eqs. ([5} that 
the linear photocurrent flows along the projection of the 
electric field onto the sample plane and it has both x and 
y components, in general. By contrast, circular photocur- 
rent flows transverse to the radiation incidence plane, 
i.e. along y axis in the chosen geometry. Despite the 
fact that the photogalvanic effects described by Eqs. ([5]) 
require the out-of-plane component of the incident radi- 
ation, they may be important for real graphene samples 
as it follows from the microscopic model, see Sec. MI Bl 

Equations ([5]) reveal that transverse and longitudinal 
photogalvanic currents vary upon change of the radia- 
tion polarization similarly to the ac Hall (photon-drag) 
photocurrents given by Eqs. ([3]). 

Thus, photogalvanic effects described by Eqs. <j5j) make 
only additional contributions to the constants A, B and 
C in phenomenological expressions ([T]), (|2|). In the case 
that the photocurrent is driven solely by photogalvanic 
effects these constants are given by, 

A oc Xc and — 2B = C cx xi- (6) 

It follows from Eqs. © and (O that the phenomeno- 
logical theory, which is based solely on symmetry argu- 
ments and does not require knowledge of the microscopic 
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processes of light-matter coupling in graphene, describes 
well the polarization dependences of the photocurrents 
presented in Fig. Q]and fitted by Eqs. ^ and ([2]). The 
incidence angle dependences presented in Fig. ^ are also 
in line with phcnomenological description. 

Hence, the phenomenological analysis is presented, 
which yields a good agreement with the experiment. 
Schematical illustration Fig. [6] as well as Eqs. (|3]) and 
([5]) show that both the ac Hall effect and photogalvanic 
effect have almost the same polarization and incidence 
angle dependences. Therefore, the analysis of polariza- 
tion and incidence angle dependencies of the photocur- 
rents is not enough to establish their microscopic origins. 
Therefore, extra arguments based on microscopic model 
arc needed. 



B. Microscopic mechanisms 

Before turning to the presentation of the microscopic 
models, let us introduce the different regimes of radiation 
interaction with electron ensemble in graphene depending 
on the photon frequency, lo, electron characteristic energy 
(Fermi energy), Ep, and its momentum relaxation rate 
1/t. We assume that the condition Epr/h^> 1 is fulfilled 
(which is the case for the samples under study) making 
possible to treat electrons in graphene as free. 

If photon energy is much smaller compared with elec- 
tron Fermi energy, ftio <C Ep, the classical regime is re- 
alized. In this case the electron motion can be described 
within the kinetic equation for the time t, momentum p 
and position r dependent distribution function f(p, r,t). 

An increase of the photon energy makes classical ap- 
proach invalid. If Hlo ^ 2Ep the direct interband tran- 
sitions are not possible and the radiation absorption as 
well as the photocurrent generation are possible via in- 
direct (Drude-like) transitions. It is worth to mention 
that if H/t -C fuj ^ Ep the transitions are intraband, 
while for Ep < hw ^ 2Ep the initial state for the optical 
transition may be in the valence band. 

In what follows wc restrict ourselves to the indirect 
intraband transitions, assuming that h/r <C fko ^ Ep 
which corresponds to our experiments with CO2 laser 
excitation. The results for the classical frequency range, 
Hlo <C Ep, relevant for THz excitation, will be also briefly 
discussed. 



1. High frequency (ac) Hall effect 

The microscopic calculation of the ac Hall effect in the 
classical frequency range, where frio <C Ep was carried 
out in Refs. PUHH. Thus, we give here only the fi- 
nal result of this work obtained within the framework of 
the Boltzmann equation with allowance for both EB [ac 
Hall effect) and qE 2 (spatial dispersion effect) contribu- 
tions. The circular photocurrent is given for degenerate 



electrons by 
3 A = A6 sin 2(p 



27rft a (H-w a 7f) d 



T9 

1 + — 



I + w 2 r| ' 



(7) 



Here q = u)/c, we have replaced for the small inci- 
dence angles q sin #0 w q9o , V is the electron velocity 
in graphene, T\ and T2 are the relaxation times of first 
and second angular harmonics of the distribution func- 
tion describing the decay of the electron momentum and 
momentum alignment [lol [fH [35j . and r = dlnri/dlne 
(e is the electron energy). The frequency dependence is 
presented by a solid curve in Fig. [3J At low frequencies 
lot <C 1 the parameter A and, correspondingly, the cir- 
cular photocurrent raises with the frequency increase as 
lot. In the high frequency regime lot 3> 1, by contrast, 
the circular photocurrent related with CacHE drops as 

1 h 

oc — , - <t:hio <s:E F . (8 

L0 6 T T 

Calculations show that for our n-type structures the con- 
stant A describing CacHE photocurrent is negative in the 
wholcf frequency range, achieves its maximum absolute 
value for lot ~ I and describes well the experiment at 
least at low frequencies (see Figured]). 

The solution of the Boltzmann equation also yields lin- 
ear photocurrents in longitudinal {jb> and jc) and trans- 
verse (js) geometries [HI EH- These photocurrents are 
proportional to the constants T\ and T2 in Eqs. ([3]). They 
describe well polarization dependences presented in Fig. [I] 
providing the polarization independent longitudinal pho- 
tocurrent as well as photocurrent contributions varying 
with the change of degree of linear polarization as sin 4ip 
and cos4<^. These constants T\ and T2 as functions of 
frequency diverge as 1/lo at lot — > and decay as 1/lo 3 
for lot 3> 1. As a result, 



I h 

3b,3col—, - < hLo < Ep. 

1 . 1^ T 



(9) 



It should be noted that the longitudinal linear photocur- 
rent can change its direction as function of the radiation 
frequency depending on the dominant scattering mecha- 
nism [la ]. 

To present a complete picture of the photocurrent for- 
mation due to Drude absorption we turn to the quan- 
tum frequency range and assume that hio ^ Ef, while 
lot ^> 1. The absorption of the electromagnetic wave 
in the case of intraband transitions should be accompa- 
nied with the electron scattering, otherwise energy and 
momentum conservation laws can not be satisfied. The 
matrix elements describing electron transition from k 
to p state with the absorption {M^ h ^ q ) and emission 

(M°™ ' q ) of a photon with the wave vector q are calcu- 
lated in the second order of perturbation theory as 



abs.q _ 
p,fc 



V p,k+q rL k+q,k 



£ k £ k+q 



TlLO 



p,p—q p—q,k 

°fc p— q 



(10a) 
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FIG. 7: Schematic illustration of the processes responsible 
for the drag effect in the quantum frequency range under 
intraband transitions (ho ^ Ef). Solid/red arrows denote 
electron-photon interaction, dashed/blue arrows denote elec- 
tron scattering caused by impurities or phonons. Filled/gray 
area shows the part of energy spectrum filled with electrons. 
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(13b) 



Here e p = e k + h-O- It is noteworthy that Eqs. (fT3"|) are 
valid provided Hlo < Ep. We note that although the 
scattering rates are not explicitly present in Eqs. (|13[) . 
the scattering processes are crucial for the photocurrent 
formation. 

If the photon frequency becomes much smaller as 
compared with the electron energies, Hlo <C £fc,£ P , but 
ujt\ , lot 2 S> 1 the photon drag effect can be described 
classically. One can check that, in agreement with 
Eqs. ©, Eqs. (HD yield 



u=± 



p,k — q k-q,k 

4 zl- q fa 



P,P+q p+q,k 
£ k ~ £ p+q 



(10b) 

Here superscript v enumerates conduction band (v = +) 
and valence band (v = — ), respectively, R k ± q k * s the 
electron-photon interaction matrix element, V p u k is the 
matrix clement describing electron scattering by an im- 
purity or a phonon. We note that the incident elec- 
tromagnetic wave is assumed to be classical, hence the 
electron-photon interaction matrix elements are the same 
for the emission and absorption processes, M™^ 1 *' 9 ~ 

^k h p 9 = ^Ipk because the number of photons in this 
wave is large. It was assumed also that the graphene is 
n doped so the initial and final states lie in the conduc- 
tion band. The intermediate state, however, can be in 
conduction or in valence bands, sec Fig. [7J 

The dc current density can be calculated as [36| 



8tt 



3 = e -f-^2[ v P T i( e p) - Vkn(e k )}\M% 



q i2 
fcl 



k.p 



[f(e k )~f(e p )}S(e p -e k -hio), (11) 



where v k is the electron velocity in the state with the 
wave vector fe, T\(e k ) is the momentum relaxation time, 
f(s k ) is the Fermi-Dirac distribution function, e k = hvk 
is the electron dispersion in graphene. 

Let us assume that the electron scattering is provided 
by the short-range impurities acting within given valley, 
intervalley scattering processes are disregarded. The ma- 
trix elements for the impurity scattering are given by 



n + = ^ 

P k 2 



I _|_ gUvie-Vp) 



V. 



pk 



Vo r 

2 



1 



v -k + - * 

pk 2 



(12) 



where Vq is real constant. As a result, one can express 
the coefficients T\ and Ti describing linear photocurrent 
in the following form (lot 3> 1) 



Ti 



-eV^L^[ /(£fc ) 



chio 4 



(£fc + £p) 2 ' 



(13a) 



CUJ° 



(14) 



where /' = df / de. In this frequency range values of T\ 
and T 2 are identical to those presented in [l6jj. Hence, 
linear photocurrents js,jc °c 1/lo 2 in this frequency 
range, see Eq. ([9]). Moreover, it can be shown that 
the circular high frequency Hall effect requires an al- 
lowance for the extra scattering and T\ oc 1/lo 4 making 
]A oc 1 /lo 3 in agreement with Eq. © . Therefore, the fre- 
quency dependence of the circular photocurrent, ja, is 
non-monotonous with the maximum at lot ~ 1. This is 
exactly the behavior observed experimentally, see Fig. [31 
where the coefficient A is plotted. Its absolute value first 
increases with the frequency and afterwards rapidly de- 
creases. Overall agreement of the experimental data in 
sample 2 (shown by the points) and theoretical calcula- 
tion (solid line) shown in Fig. [3]is good. The theory, how- 
ever, does not describe the abrupt frequency dependence 
and change of the photocurrent's sign observed in sam- 
ple 1 (see gray circles in Fig. [3]) . In order to understand 
this behaviour we analyze the possible contributions of 
photogalvanic effects. 



2. Microscopic mechanisms of photogalvanic effects 

Real graphene samples arc deposited on substrates. As 
we already noted above, it results in a lack of an inver- 
sion center and, correspondingly, allows for the photogal- 
vanic effects. Phenomenological analysis demonstrated 
that the polarization and incidence angle dependences of 
the photogalvanic current are almost the same as for the 
oc Hall effect. It follows from the general arguments and 
phenomenological considerations summarized in Eqs. ([3]), 
that the photocurrent can be generated only with al- 
lowance for z-component of the incident electric field. 
However, for strictly two-dimensional model where only 
7r-orbitals of carbon atoms are taken into account, no 
response at E z is possible. Therefore, microscopic mech- 
anisms of the photogalvanic effects in graphene involve 
other bands in electron energy spectrum formed from the 
(T-orbitals of carbon atoms. 
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71— orbitals, 
PT band 




a-orbitals, 
Pt band 



FIG. 8: Schematic illustration of indirect intraband transi- 
tions with intermediate states in band which interfere 
with Drude transitions (shown in Fig. [TJl and give rise to the 
photogalvanic effect. Solid red arrows are electron-photon in- 
teraction, dashed blue arrows are the electron scattering. 



There are 6 irreducible representations P } , P 1 ~ , P%, 
P 2 ~ , P 3 , and P 3 ~ at K (or K') point of the graphene's 
Brillouin zone. The conduction and valence band states 
transform according to the P 3 ~ representation: there are 

two basis functions p^ z \ p z 2 ^ being odd at the reflection 
in the graphene plane z = 0. Symmetry analysis [37l l38j 
shows that the transitions in z polarization arc possi- 
ble between these states (transforming according to P 3 _ ) 
and the states transforming according to P 3 + . The lat- 
ter representation is described by two functions and 
s^ 2 ) which do not change their signs at the mirror reflec- 
tion z — > — z. Under the symmetry operations which do 
not involve z — > —z these wave functions transform like 
Pz^ , Pz ■ Representation P 3 + corresponds to a orbitals of 
carbon atoms which form remote valence and conduction 
bands of graphene. Microscopic calculations performed 
within the basis of 2s and 2p atomic orbitals j37H39| show 



that the distance from the P, 



states forming conduction 

+ 



and valence bands and closest deep valence bands P 3 
A, is about 10 eV. It is remarkable, that the electron 
dispersion in these bands has the form, similar to that 
of conduction and valence bands: i.e. energy spectrum 
near K (or K') point is linear, however, with different 
velocity, as it is schematically illustrated in Fig. [5] 

Microscopically, circular photogalvanic effect arises 
due to the quantum interference of the Drude transitions 
represented in Fig. [7] (for q — 0) and the indirect intra- 
band transitions with intermediate states in P 3 + bands 
depicted in Fig. [5J similarly to the orbital mechanisms 
of the photogalvanic effects in conventional semiconduc- 
tor nanostructures (40M421 ]. Indeed, matrix elements of 
Drude transitions are proportional to the in-plane com- 
ponents of electric field E\\ and electron in-plane wave 
vectors in the initial k, and final p states. The matrix 
elements of the indirect transitions via P+ band are pro- 
portional to E z and do not contain linear in fc, p contri- 
butions. As a result, the interference contribution to the 



transition rate is proportional to both E\\ and E z and 
to the in-plane wave vector components giving rise to dc 
current. The presence of the substrate allows electron 
scattering between the states transforming according to 
P 3 + and P 3 ~ representations: for instance, the impuri- 
ties located near the substrate surface or the phonons, 
propagating in the substrate, or the impurities adsorbed 
from the air to the graphene create an effective potential 
which is not symmetric with respect to z — > — z mirror 
reflection. Hence, the interference contribution to the 
transition rate is non- vanishing. 

Let us denote a orbital states transforming according 
to P 3 + orbitals as +' and — ' (we recall that the super- 
scripts + and — denote the conduction and valence band 
states in Eq. (fTU|) . respectively). We assume that the rel- 
evant interband optical matrix element has a form p3j 



P+ + — _ R++ 
n kk ~ n kk 



= A z ip , 



(15) 



where ipo is the momentum matrix element between a 
and 7T orbitals, po is assumed to be real (and the momen- 
tum matrix element is imaginary). 

We also need to define the form of the interband scat- 
tering matrix elements. We have already noted that the 
phonons in the substrate or the impurities positioned 
either above or below the graphene sheet can provide 
the scattering between the bands transforming by P 3 + 
and P 3 ~~ representations. In addition, the impurities or 
phonons should also provide the scattering within the n- 
orbital band. Such a scattering should be short-range in 
order to allow the electron transition between <j and tt 
orbitals. We assume that the interband scattering also 
takes place between the similar combinations of the Bloch 
functions. We take the scattering matrix elements in the 
following form for the interband scattering for the rele- 
vant processes [43j : 



v p V = v p V 



Vi 



1 _|_ e i(i/>fc-ip P ) 



(16) 



with Vi being the real constant. 

The second-order matrix element for the scattering- 
assisted optical transition via a orbital can be written 
as 



V++ p+ -t 

v pk ^kk 



R++'V + , + 

PP pk 



■k,+ 



■ k ,+' + huj e k ,+ - e P ,+> 



(17) 



Here Sk, v with u = + or +' describes electron dispersion 
in a given band. Corresponding processes arc depicted in 
Fig. [5J To simplify the calculations we assume that the 
dispersions of electron in a and ir bands are the same. 
The allowance for difference of effective velocities will re- 
sult in the modification of the results by the factor ~ 2. 
Equation ([T7| under assumption that A ^> fiio, Ep trans- 
forms to 



M, 



pk 



.eA z poV! 

1 

2moc 



2hw 
A 5 "" 



(18) 
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It is the quantum interference of the transitions via a 
orbitals described by Eq. ([15]) and Drude transitions de- 
scribed by Eq. ([TO]) (where one has to put q = 0) [13, EH 
that gives rise to the photocurrent. The photocurrent 
density under the steady-state illumination can be writ- 
ten as [cf. Equation ([II]) and Ref. ^} 

3 = ^-Y. m {K^K^ v ^{e v )-v k r 1 {e k )\^ 



k.p 



{f(e k ) - f(e p )]S(e p - e k - fej). (19) 



Making necessary transformations we arrive at the fol- 
lowing expression for the constant \c describing circular 
photogalvanic effect: 

4?rw \ ^ Ti{e p )e k + Ti(e k )e p ^ 
Xc — — 7— / ; x 

lf(e k ) - f{e P )]5{e p -e k - foj), (20) 



where 



2ire 2 vp (Wi) 



JJIqCLJ A z 

and (. . .) denote the averaging over disorder realizations. 
Equation ([20| is valid provided lot ^> 1 and hw < Ep. 
The treatment of the general case is given in Appendix 
to the paper. 

The direction of the current is determined by the sign 
of the product (VqVi) and the radiation hclicity. The 
averaged product (VoVi) has different signs for the same 
impurities, but positioned on top or bottom of graphene 
sheet. It is clearly seen that the photogalvanic current 
vanishes in symmetric graphene-based structures where 
(V o V 1 ) = 0. ' 

In the case of the degenerate electron gas with the 
Fermi energy Ep and in the limit of huj <C Ep Eq. (|2"0)) 
can be recast as 



Xc 



aed (V V 1 )E I 



A {V a 2 } Hu>' 
we introduced effective dipole of interband transition 

ep h 



(21) 



mo A 



In Eq. ([21]) a is the fine structure constant. It follows 
from Eq. ([2"Tj) that the circular photocurrent caused by 
the photogalvanic effect behaves as 1/cj at lot ^> 1, 
Hlo <C Ep, i.e. it is parametrically larger than the cir- 
cular ac Hall effect which behaves as 1/w 3 , see Eq. ([8]). 
This important properly is related with the time reversal 
symmetry: the coefficient \c describing photogalvanic ef- 
fect is even at time reversal while T\ describing caHE is 
odd. Therefore, circular photocurrent formation due to 
photogalvanic effect is possible at the moment of photo- 
generation of carriers, making extra relaxation processes 
unnecessary 



As discussed above experimental proof for the CPGE 
comes from spectral sign inversion of the total photocur- 
rent observed in sample 1 [see Figs. [3l [4] and EJa)]. Let 
us estimate the circular photocurrent and compare it to 
experiment assuming that the photocurrent in sample 1 
is dominated by the CPGE. Taking d = 1 A, A = 10 eV 
we obtain 

, (VoVi) E F 1n _n A cm h . 

(22) 

In the studied frequency range of CO2 laser operation 
Ep /(Too) 3. Considering the strongly asymmetric scat- 
tering, where (VoV!)/(Vq 2 ) ~ 0.5, our estimation yields 
A ~ 2 x 10 -11 (A cm)/W which is in a good agreement 
with experiment [see Figs. [3]. The values of the circu- 
lar photocurrent driven by the PGE and by CaCHE are 
similar for Tilo <~ 100 meV. It means, that for lower fre- 
quencies, the CaCHE dominates, since it has stronger 
frequency dependence, while for higher frequencies, the 
circular photogalvanic effect may take over. While the 
sign of the circular ac Hall effect is determined solely by 
the conductivity type in the sample and the radiation 
helicity, the circular photogalvanic current sign depends 
on the type of the sample asymmetry. In general, these 
two effects may have opposite signs which may result in 
the sign inversion observed in experiment, Fig. [3J 

The strongly asymmetric scattering might be exactly 
the case for the short range impurities positioned on the 
substrate surface or adsorbed from the air on the open 
surface of the sample and which provide the same effi- 
ciency of both inter- and intra-band scattering. Obvi- 
ously, the degree of asymmetry and even its sign, which 
reflects the coupling of the graphene layer with the sub- 
state, depend on the growth conditions and may vary 
from sample to sample. This explains the fact that the 
sign inversion is detected only in some studied samples. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have carried out the detailed exper- 
imental investigation of the photocurrents in graphene in 
the long wavelength infrared range. The photocurrents 
were excited by pulsed CO2 laser at oblique incidence in 
large area epitaxial graphene samples. The magnitudes 
and directions of the photocurrents depend on the ra- 
diation polarization state and, in particular, the major 
contribution to the photocurrent changes its sign upon 
the reversal of the radiation helicity. 

Phcnomcnological and microscopic theory developed in 
this work show that there are two classes of effects being 
responsible for the dc current generation driven by po- 
larization of the radiation. Firstly, the photocurrent may 
arise due to the joint action of the electric and magnetic 
fields of the electromagnetic wave (or transfer of the ra- 
diation wave vector to the electron ensemble). Secondly, 
the current may be generated due to the photogalvanic 
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effects which become possible when the inversion sym- 
metry is broken by the presence of the substrate. In this 
case, the magnetic field of the radiation or its wave vector 
are not important, but the asymmetry of the structure is 
needed. Arguments based on the symmetry to the time 
reversal show that even in the case of small asymmetry 
of the sample, the circular photogalvanic effect can be- 
come paramctrically dominant at high frequency due to 
weaker decrease with an increase of the frequency (1/oj 
as compared with 1 /w 3 for CacHE) . While both types of 
photocurrents are indistinguishible on the phenomeno- 
logical level, investigation of their frequency dependence 
allowed to distinguish them and provided direct exper- 
imental proof for the existence of CPGE in graphene. 
Microscopic theory of ac Hall effect and CPGE give a 
good qualitative as well as quantitative agreement of the 
experiment. 

Our experiments also demonstrated that photocurrent 
exhibits resonance behaviour at frequency close to edge 
of the reststrahlen band of the SiC substrate at about 
huj rj 121 meV. The resonance is observed for all pho- 
tocurrent contributions and may indicate an importance 
of the graphene coupling to the substrate and role of the 
phonons in the substrate. The origin of the resonance 
remains unclear and determination is a task of future 
work. 



in asymmetric structures z component of the incident 
electric field gives rise to the temporal oscillations of the 
electron momentum scattering time T\{t). As a result, dc 
current is formed 



j OCT! (t)E\\(t), 

where over line denotes temporal averaging. 

The method developed in Ref. [42j can be generalized 
for graphene. Indeed, the processes depicted in Fig.[8]and 
described by the matrix element (|18p can be interpreted 
as the E z induced correction to the electron scattering. 
Equation (fT8|) can be recast as: 



eE z (t)PaVi 
2mo 



1 _j_ gi(Vfc-Vp) 



As a result, the correction to the electron momentum 
scattering rate is given by 



i i i 



27T 



pk 



S(Sp - £fe)[l - COS (ip p - Cfk)] = 

CeE,(t), (A2) 



Acknowledgments 

We thank E. L. Ivchenko, S. A. Tarascnko, 
V. V. Bcl'kov, D. Weiss and J. Eroms for fruitful dis- 
cussions and support. Support from DFG (SPP 1459 
and GRK 1570), DA AD, Linkage Grant of IB of BMBF 
at DLR, Applications Center "Miniaturised Sensorics", 
Swedish Research Council, SSF, RFBR, Russian Min- 
istry of Education and Sciences and "Dynasty" Founda- 
tion ICFPM is acknowledged. 



Appendix A: Photogalvanic effects in classical 
frequency range 

In the case of Huj <C Ep photogalvanic effects allow 
a simple and physically transparent interpretation [42j: 



where 



C = s 



(VoVi) d e k 
v 2 h A : 



(A3) 



where S is the sample area. 

Following Ref. we obtain the photocurrent density 
in the following form: 



8aed £F {V Vi) 
3 = — r \„ n / Ix 



hA (VJ) 
e\\e* + ete 



1 + (ujry 



+ i(eue z - ete z ) 



UJT 



(A4) 



In agreement with symmetry considerations, Eq. (|5|), 
both linear and circular photocurrents are allowed. For 
ujt > 1 Eq. (|A"4|) agrees with Eq. ([2i~]) . 
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